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Abstract 



The potential between sources in arbitrary representations of the gauge group is studied on an anisotropic 
lattice in a spherical model approximation. It is shown analytically that for half-integer j and j' in the 
confinement phase the potential rises linearly , whereas for integer j and half-integer j' it rises infinitely 
which means a strong suppression of the combination of such states . For integer j and j' the potential 

fvj ■ shows Debay screening and Coulomb behavior in the deconfinement phase . It is also shown, that ( x ) ^ 

'^ ■ (x) "* when (x) -^ 1 and is in agreement with the mean field theory prediction, and (x ) ^ ix) for 



m . 

("■^ , (x) ^ 1 which agrees with MC experiment. String tension model-computed for sources invariant under 

("^ ' center group transformations demonstrates Casimir scaling in the intermediate distance regime and turns 

\^ ' into zero at large distances. 
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1 Introduction. 

Statistical QCD evaluated on the lattice by means of computer simulation is perhaps the only 
case in statistical physics where the critical behavior can be calculated from the first principles 
of dynamics ||l|, g, |3|. The power of the approach is best seen in the thermodynamics of pure 
SU{N) gauge theory, i.e., for the systems consisting of gluons only [^ ^, g. The main features 
considered are the deconfinement transition and the properties of the hot gluon plasma f^ . Lattice 
gauge theory (LGT) provides a powerful tool to compute nonperturbative properties of strong 
interactions, such as a potential between pairs of infinitely heavy sources. 

Its known from Monte Carlo simulations that probes in representations insensitive to the 
centre gauge group transformation Z{N) yield a screened potential, while those sensitive to this 
transformations provide a linearly rising one F « aR (where a - string tension, R - distance 
between sources) H, H, H. As indicated in mn, the potential between sources transforming as 
arbitrary irreducible representations of the gauge group provides a good laboratory for testing the 
dynamics behind screening and confining mechanisms. In particular, one can measure the relevant 
distance scales, i.e. those that separate "short distances" (essentially perturbative physics), from 
the "long distances" . The dynamics of the intermediate region has been found quite rich in many 
aspects. Indeed, some very simple numerical simulations have shown that adjoint sources, for 
example, may feel a linearly rising potential at intermediate distances |ll], |l2|, |l3[ |lj, |l^. Such 
sources appear to "deconfine" at precisely the same critical temperature Tc at which fundamental 
sources deconfine [^^ |l^. The string tension is representation-dependent and appears to be 
roughly proportional to the eigenvalue of the quadratic Casimir operator pil O, nsl O. 

As the distance between color sources with zero " n-ality" increases they are eventually screened 
by gluons. This asymptotic regime extends from the color-screening length to infinity, and in the 
case of SU{2) gauge group the string tensions become equal for half-integer j -s and turn into 
zero for integer ones. In particular, the string between quarks in an adjoint representation must 



break at some distance that presumably depends on the mass of "gluelumps" (i.e. on the energy 
of a gluon bound to a massive adjoint quark) p9| |. 

^From universahty hypothesis it follows, in particular, that all higher representations sensitive 
to the center of the gauge group should be equivalent order parameters with the same critical 
behavior as the fundamental representation. However different behavior was found for the higher 



representations at SU{2) lattice gauge theory in (3 + 1) dimensions |l^, 17|. 

The purpose of this work is an analytical investigation of the potential between sources in 
arbitrary representations of the gauge group near the critical point on an anisotropic lattice in a 
spherical model approximation. 

The application of spherical model in statistical physics has a long history since the time it 
had been introduced to investigate critical phenomena in the ferromagnet [EQ| and until recently 
(see, e.g., [^ ^)- Stanley ||2^ established the correspondence between the spherical and the 
Heisenberg models. 

Although this model is of no direct experimental relevance, it may provide useful insight since 
many physical quantities of interest can be exactly evaluated with its help. In this context, the 
spherical model is quite a useful tool in providing explicit checks of general concepts in critical 
phenomena, see [|4| |2^, |2^, g^l ■ Recently ^ it was successfully used for studying the transitions 
between a paramagnetic, a ferromagnetic and an ordered incommensurate phase (Lifshitz point). 
Models of this kind were investigated extensively (see recent review in |28| ). The spherical model 
predicts reasonable values for critical exponents. Moreover, a 'basic' set of exponent relations is 
also satisfied by spherical model for 2 < d < 5 p9| |. 

2 LGT on anisotropic lattice. 

The simplest way to increase the accuracy in numerical study of LGT at large distances and high 
temperatures is to use anisotropic lattices in which the temporal spacing is much smaller than that 



in the spatial directions (^ ~ a^/ar ^ 1, the Hamihonian hmit). This approach becomes a very 
popular technique at the present time since it allows Monte Carlo calculations to be carried out 
with reasonable computational resources and reap many of the benefits of fine lattices, while still 
using cheap coarse spatial ones. There are convincing arguments pO| , p31| that in the anisotropic 
Wilson gauge action no coefhcients are to be tuned in order to restore space-time exchange 
symmetry up to O(a^) errors. On such lattices heavy quarks will not suffer large lattice artifacts 
as long as their masses are small in the units of a-r- 

To compute the partition function on a lattice of a size Nr x N^ {Nr is the temporal extent, 
N^ is the spatial extent of a lattice) after |Q (see also [^ Q ) we use the anisotropic lattice 

Sg — PtSe + PaSM, (1) 

where Se is electric and Sm magnetic part of the gluodynamics action and 

~ ~ 2N 

f3r = K (C, P) » Pa = Pli (e, /?) , P=—. (2) 

with 



i = aajar. (3) 

In the weak coupling region Pt,<jI P ~ 1 + O (1//9) |3^i HI and in this case ^ does not differ 
essentially from ^. Moreover, in the "naive" limit 

U, (x) ~ 1 + la^g.Al {x) T, (4) 

and one must put exactly ^ = ^. However, in LGT the effective values of the integration variables 
Uy [x) are far enough from the "naive" limit (especially in the deconfinement region). The devia- 
tion of ^ from ^ is regarded as 'quantum corrections' or 'renormalisation'. The behavior of ^ (^, /?) 
in the vicinity ^ ^ 1 has been studied both perturbatively p7| and by non-perturbative methods 
|3S| , pO| , PI , ^ |4^ . Since we cannot exclude, that in the area of the asymptotically large ^ the 



behavior of ^ may be essentially different, we will not specify the dependence of ^ on ^, but only 
assume that ^ (^, /3) may be made arbitrary large by increasing ^ at fixed (3. 

It is well known 1 44, 451 (see a review in p6[) that there is a variety of ways to get (in accordance 



with a universality hypothesis |32[ ) the effective action 

3 Nl 

-5e//=7(/9,C,A^r)Re^^X.x:+^ (5) 

1=1 X 

expressed in terms of the Polyakov lines Xx — Tr J^^Jj^ C/q {x,t) in j— representation of SU{N) 
group (we omit index j for fundamental representation). The crucial assumption for obtaining 
(||)is to discard the magnetic part {Sm ^ 1/0 of action. Then Sq [U] ~ Se with 

3 

Se ^-P£,Y.Y1 ^^ Tr {c/o {x) U, {x + 0) ul {x + i/) C/t (x) - l} . (6) 

X I/— 1 

As a rule it can be done in a strong coupling approximation [E3, Bs]. However, as it is 
pointed out in [p2| at high 'temperatures' Tsy = £,/ (arNr) = £,T and for the couplings one gets 
Pa- ^ l/TgY a-nd (3^ ^ Tsy. Hence we may assume that for finite T the magnetic part of action 
may be neglected |33] compared to the electric one {Se ^ i^ a, Hamiltonian limit (^ ^ 1) even 
at low temperatures T. Such an assumption doesn't look as harmless as in a strong coupling 
case. Indeed, there are serious reasons to believe that the magnetic part of the action may be of 



crucial importance for creating the confining forces |^, 48, ^, even at high temperatures [pO|| . 
Therefore, strictly speaking, QCD without a magnetic part may be considered in a weak coupling 
area only as a specific ('toy') model. 

Summarizing for the partition function with static sources i^x we may eventually write 



Z~y exp^Re 7 ^ XxX:+, + ^ »7x • X^ •^H'^^- (7) 



3 SU{2) gauge model 



In a simple case of SU{2) gauge group one may write 



d^. = Bin (- j — 



27T 

; X^ = Xx = 2cos— , (8) 

— 27r 



so by introducing new variables: cTx — 2cos^ and ax = 2sin^ the partition function (^ can 
be rewritten as 

Z^ exp < 7^o-^crj;+; +^77^0-2; > Jj5(cr^ +5^ -4)dcrj;fT^dfTj;. (9) 

"'-°° [ X,l X j X 

To compute the partition function we use the spherical model approximation^ : 

(10) 



X \ X / 



c+ioo 7 
c—ioo ^""^ 



where constant c is chosen so that the integration path is placed at the right side of all singularities 
of the integrand to ensure the legitimacy of interchanging the integration order over ds and daxdcTx- 
Therefore we can rewrite (^ as 

Z = const I dse^"^ I exp {-axAxx'Tx' + VxCTx) d(7x e~'"^''aldax, (11) 

J J —00 J —00 

where 

d 
Axx' = sSx,x' ~ l^^Sx+l,x'- (12) 

1=1 
It is clear that integration over ax can be done trivially 

f'[le-'~^^-ald'Sx ^eyii>l~N^^ln{s-3-f) + const\. (13) 

Integration over ax may be fulfilled after Fourier transformationQ 

ax=Y.Cke''"',qi = ^; / = 1;2;3; k ^ 0, ...,N, - 1 (14) 

k 



1 We want to stress, that the spherical model approximation does not break the global Z (2) - invariance of the 
model. 

^ Discrete variables qi for N^ S> lean be considered as continuous: < g; < 2it. 



which diagonahzcs bihncar form 

Y.<J.A^.'CJ^'=NlY,\Ck?A{q). (15) 

x^x' k 

where 

d 

A{q)^ s-j^cosqi. (16) 

1=1 

Now we can do the integration over (k and write for partition function 

Z = const / dscxp Y^ VxA-^,r]^, + TV^ Us - 6 (s) - - Ins j , (17) 



where 



b{s) = isp{lnA} = l£^ (S)'^^ fs -7X^cos(Z/ j , (18) 



and 



3 

r27r 6^*93(^-1') f ia' 



37^.-,'-/ 1^3^ ' --R3^-;x-J, (19) 



37 3 



with 



'''i^^^ri '"'*'H5)^nii*' <'"* 



and /„ (x) - modified Bessel function. 
Taking into account that 



x y 

we obtain for the uniform sources rjx = ??o 



2^ ^2:y - 2^ ^a:y ^ ^ _ 3^ ' (^1) 



with 



Z ~ consi / ds exp {iV$ (s)} , (22) 



<J>(s).4s-Hs)-|lns+^^. (23) 



To calculate ( p^ ) we use the steepest descent method that in the limit N^ -^ oo gives an 
accuracy of order O i^j. The saddle point sq is determined by $' (sq) = and considering that 



one may write |g9|] for |2. — 1 = e ~ 0. 



we get 



B3(J + l;0)»|-^\/|^v^, (24) 



■^'-*-m-4M-.;r7^-TTA^^ m 



2 I 27 ^7 / 27(l+e) 4(87) e^ 



Hence for 7 < 7c the solution of (g5|) for ?/ = is given by 

so ^37(1 + ^(8717)2(^^-1^ ^(7c-7)J, (26) 

where the critical coupling -fc is defined as the solution of $' = at the lowest value of sq -^ 87 + 
(and 77 -^ 0), which evidently gives 7c ~ 3/16. In a vicinity of the critical point one may write for 

It is easy to see, that in the deconfinement region (7 > 7c) the solution of $' = with 77 = 

is absent. The only way to regain the solution in such region is to tend sq — )■ 87, because in this 

2 
case the term -. — ^° .^j must be preserved even in a limit rjo ^ 0. Now the solution 

so ~ 37 + -fllL= -^ 37 (28) 

/l _ 2£ 

7 

sticks to the point 87 so in this region we can write 






and 



(x) 






^2h-^) 0(7-7c). (30) 

770^0 



Hence in the deconfinement region (x) ^ but steadily approaches zero when 7 ^ 7^ + in 
accordance with the weU known fact, that SU{2) - gluodynamics undergoes second order phase 
transition. 

For the two point correlation function for i? 3> 1 we get 

(X0X.> = ^ + (x>^-c^+4(l-^).(l-^), (31) 

where c = const and string tension (in lattice units) is given by 

a = 3,(l-i)«(l-i). (32) 

The expression for correlation function ds^) as (EQ) is just the Fourier transform of the 
propagator \D{q) + ?7i^] where D{q) = X]i=i 4sin^ ^ is the lattice laplacian in d-dimensional 
momentum space. In the perturbation theory expression (GO) can be attributed to the lattice 
one-particle exchange. A basic non-perturbative feature of the high temperature plasma phase of 
QCD is the occurrence of the gluon screening mass ||5l|], so parameter a, defined by ([3^ ) can be 
interpreted as screening mass as well. The expression (|20| ) is not specific for the spherical model 
and appears in many popular approaches, for example in the matrix model (see e.g., |52| ) and at 
the perturbation expansion (g^ ^ 1) (see [|3[ p^). 

A more detailed case is considered in the Appendix. Here we present only the final results : 
We work in the parameter area 7 ~ 7c and R ^ 1 where^lQ Jj /^oo "^ ^'^ ^'^'i (x) ^ 1 so we 
preserve only the terms with lower powers of Aq ^/Aq^q and (x). 

It is easy to see that in all cases the correlation function 

C(^'(^"')^(x^^\i^"^)-(x(^"))(x(^"')) (33) 



J &/ 


(x^^^r^) 


both integer 


(x(^-) ) (x(^') ) exp { ixf A-],c, + if J i) (a- ) " c, } 


both half- int. 


(x(^))(x(^')) + 'z(/;j)^ 


j int.; j'half- int. 


(x(^-))(x(^')) + (x(-')>(x(^'-^>A->3 



Tab. I Pair correlation for states in higher representations 

exponentially decreases in the confinement phase ((x) = 0) and can be used for the computation 
of a, which plays a role of string tension when (x ) (x^"' ) = ^^^ may be interpreted as a 
screening mass for (x'"''') (x^"* -^ ) 7^ 0- 

4 Conclusions. 

As it can be seen from TABLE 1, our model gives linearly rising free energy F ~ aR in the 
confinement phase for half-integer j and j' . The potential between states with j integer and j' 
half-integer (or vice versa) turns into infinity which means strong suppression of such combination 



of states. For integer j and j' we get the screening potential: F • 



In the dcconfincmcnt 



phase ((x) ^ 0) we get F ^ —l/R. We would like to stress here that in the deconfinement phase 
((x) 7^ 0) for integer j and j' the screening mass is twice as much than in confinement ((x) = 0) 
one. 

As it was predicted by the mean field theory [p5l 



(3) 



(x) 



2j 



(34) 



This result was confirmed in MC simulations ||56| and as it can be seen from TABLE 1 roughly 
agrees with such data when (x) }t 1. However, as it was pointed out in [ p7| in the vicinity of 
critical point the relation (^^ does not agree with MC experiment. TABLE also shows, that our 
result (Egl) 



U) 



(x), 



(35) 



obtained in such area agrees with the data given in p7[ and with theoretical predictions developed 
both in framework flux picture in a strong coupling approximation and within field theory ((/)*) 
on the basis of universality arguments p7| . 

Finally, there is the issue of the Casimir scaling (whose importance were emphasized in 
[|l4 [l2| ) of higher-representation string tensions in the intermediate distance regime. For val- 
ues R far from asymptotic region {aR < 1) we can write A^]^ « cq • (1 — aR) and (x'""'''x^"' -^ ) — 
(x ) ( X^"* ) 6xp |— a • (j'j — j] cR + const\ where c — 2c2Cq is independent either from R or 
from j and j' so it follows from ( |54|) and (|5q): 

a(j;j')^«-(/j-^)- (36) 

Therefore, the adjoint sources may 'feel' a linearly rising potential at intermediate distances, 
and appear to "deconfine" at precisely the same critical point at which the fundamental sources 
really deconfine. This qualitatively agrees with Casimir scaling observed in such parameter region 
[|ll| Q . It is easy to see, that at very large distance scales Casimir scaling breaks down and color 
screening is set. 

As it was pointed long ago, sources trivially transformed under centre gauge group transfor- 
mations can not be used as an order parameter. An average value of a pair of such sources also 
is not a good laboratory for critical phenomena studies, because the magnetization Mv-''^ -' = 
lim/v^oo ( Xo Xr ) = (x'^"''') (x^"' w differs from zero and masks the Debay term. Therefore 
^''r ~ \ Xo Xr ) ~ (x ) (x^"' -^ ) ^ with i? ^ cx3 so Xo and xk ^^^ statistically indepen- 
dent at asymptotically large distances and long range order is absent either for (x*-"'-') (x^"' -^ ) = 
or for (x'-"''') (x^"* ' ) t^ 0- On the other hand, the fact that the correlation function Cj^-' exponen- 
tially decreases for all j and j' can be used for determination of string tension even in the cases 
when (x ) or ( x^"* ' ) arc invariant under Z{N) transformations. 

We also wish to note that computing the Binder cumulant B, defined as [p8[ - 

10 



we get for a case SU{2) gauge group 

Therefore, in a vicinity of a critical point we get i? « because (x) = in the confinement region 
and in deconfinenient region, if we are close enough to the critical point. Vanishing of B means 
that the distribution of x is Gaussian, therefore the spherical model must give reliable results in 
both regions. 

5 Appendix 

The n - order correlations (nfe=i Xx^ ) may be easily computed after Fouirier transformation of 

Vr{£j^J"^ ur (z) e*- 1^ ; ur{z)^ p \Jr[^) c'^'^^da (39) 

and we get in a spherical model approximation 

\fc=l / -^""^ \ { X ) I k=l k=l 

pc+ioc |.oo n ^ 

^ ds e-npi-ffxAx.x'O-x' +r]xCrx}Y\_'^2j^{zk)-;T—Y\_d'Crx (40) 

J c—ioo J —oo 7,_i ■^'^ 

with 

n 

Vx = VQ + iYl ^fe*^^' ; Vo^O- (41) 

fc=i 

Therefore after integration over a^ 

' " \ ^^ poo " J- -1 "1 " W 



\fc=l / V v^^v.i .,-00 j,^-^ ^ J f.^^ 



and taking into account 
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770— *u 4 4 



( -^VxA^^x, 



^=«oJ^„^0 



r^fe"fc,r^'-+*(X>X!^'=' 



we obtain inverting Fourier transformation 



nx^0^nu.4f), 



(43) 



\k=l 



fc=l 



with 



Q = ^- / Q CXp <! - ^ (cTfe - (x)) ttfe.r (cTr - (x)) > J| rfcTfe , 

^' -^-^^ [ fc,r = l J fc=l 



(44) 



where ak,r is n x n matrix is defined by the condition a^. j. ~ A,^^ ^^ (■so)- In particular one can 
write 



(x") = (-i^/2)^ff„((x>/V^2^ 

^l + fc(x)V^)(^/2)*n!/(f)! [f] 
(X)(i?/2)^n!/(ii^)! [f]. 

{xr + {xr-'Dn{n-l)/4 

where the dispersion I? = ( (x — (x))" ) is given by 



(45) 



7t-l . 

2 ' 



7 < 27c, 
7 > 27c, 



D = (x') {xf = ^0,0/2 = 



1/ 


47; 


7^ 7c, 


<1 


4- 


1 . 
27' 


|7c < 7 ^ 7c, 





(46) 



The average value of the character x in ^in arbitrary irreducible representation j may be 
computed as^ 



Xo) 



m 



E (-1)' 



2j — ra 



\ zj — zm / \ 



(47) 



•^ Here [j] refers to an integer part of j and H,i (x) are Hermitian polinomials of order n. 
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or 



Aj) 



Ij] 



= E (-ir 



2 j — m 



ix''-'"-) 



(48) 



that for 7 < 27c, meaning (x) ^ 2Z) immediately gives 



ij) 



3 (2j)! 
4 






O-iY- V^ (^ 87 



(49) 



In the 'deep' deconfinement area 7 > 27c or I (x) > 2D 



M) 



„,(<!>).«. ._^ 



(50) 



The pair correlation 



Xl X2 ; - U2,, (yj U2,, (y 



vdetf 



>2ji 



o-i 



Us,,, ^ u,,, ^ 



2jl 



0-2 



(51) 



exp{- {(Jk - (x))aA;,r (ct,. - (x))}d^cr, 



where matrix a^.^ is given by 



ai,i ai.2 



V 



0^2, 1 ^2,2 



/ /l-l /1-1 \ 



("^0,0 






yl 



0,0 



A, 



0,R 



\ ^0,R ^0,0 J 



(52) 



one may do computation expanding over powers Aq ^ and taking into account simple relation 
XX — X + X • III particular for magnetization Mjj/ we get 



M,,'.(x^^\2')) = (x(^))(x(^')), 



(53) 



In the confinement area (7 < 7c) for integer ji^2 = [ji,2] for aR ^ 1 



Xl X2 



Ol)\ /^0'2) 



X" '/exp 



nj2 - 1 



2 A 



0,fl 



(54) 
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In particular for ji = J2 = j we obtain for aR ^ 1 



T'-"2 



> ~ exp < 



aRC2 

27 



(55) 



where C2 = {j + ^) ^ i is quadratic Casimir operator. 
For half-integer ji^2 == b'1,2] + 5 



(^0,0) 



(56) 



with 



and 



9(ji;j2) 



((x(-+- 



)\ + /v(^-^)\U/vO-^ + ^)\ + /v(^-^) 



+ (X 



(57) 






J-2 



(58) 



Therefore taking into account 



(x(^)x) = (x(^+^)) + (x(^-^)) - (x> ^^'^^ + (x> ^^■-\ 



(59) 



we obtain 



Ah) Ail) 

A,X\ A.X2 



U.)\ /..(.,)\ , 2(xX^^-))(x(^-)x) ,,i 



(^o-,S) 



o,it: 



(X) ^^■^+^-''^ exp<|87^((x) + ^) Aq-^^. 



(60) 



In particular we get 



(x(^Vi = (x^) 



VxA^^Q 



2^0,0 — (x) + 2^0,0' 



nii:=vo^o 



(61) 
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and 



2 a2 



(xIxr) 



_(f_d_ 
dri\ dril ' 



,\vA-'-n 



(62) 



•q^^m)^^ 



ixY 



+ ^^o:o) +^(^;.,o)' + ^^o,h(x)% 



Xi),ii) 



therefore, correlation function for characters in the adjoint representation (XoXii) may be 
written as 



xi'\^') - (xlxl) - 2 ix') + 1 = {x^'^y + \ (a-^,)\ lA,^n U? 



(63) 
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